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ABSTRACT 

An order of magnitude more dwarf galaxies are expected to inhabit the Local Group, based on cur- 
rently accepted galaxy formation models, than have been observed. This discrepancy has been noted 
in environments ranging from the field to rich clusters. However, no complete census of dwarf galax- 
ies exists in any environment. The discovery of the smallest and faintest dwarfs is hampered by the 
limitations in detecting such faint and low surface brightness galaxies. An even greater difficulty is 
establishing distances to or group/cluster membership for such faint galaxies. The M81 Group provides 
an almost unique opportunity for establishing membership for galaxies in a low density region complete 
to magnitudes as faint as Mr' = —10. With a distance modulus of 27.8, the tip of the red giant branch 
just resolves in ground-based surveys. We have surveyed a 65 square degree region around M81 with 
CFHT/MegaCam. From these images we have detected 22 new dwarf galaxy candidates. Photometric, 
morphological, and structural properties are presented for the candidates. The group luminosity function 
has a faint end slope characterized by the parameter a — —1.27 ± 0.06. We discuss implications of this 
dwarf galaxy population for cosmological models. 



Subject headings: galaxy groups: individual (M81 
galaxies: photometry - galaxies: 
colors, radii) 

1. INTRODUCTION 

Standard ACDM, hierarchical structure formation mod- 
els predict an order of magnitude more low mass ha- 
los in group environments and as satellites around Milky 
Way-like galaxies than there are obse rved dwarf galaxies 
fKlypi n et al.lllQOOl [Moore et al.|1999| . The steeper slope 
of the predicted mass function over the observed luminos- 
ity function (LF) is the source of the 'missing galaxy' prob- 
lem. There is growing evidence that there are fewer dwarfs 
than predicted in all environments, but the situation is 
even worse in field and group environments. Simulations 
predict low density environments to have fractionally more 
dwarfs than the high density cluster region s, whereas ob- 
servationally the reverse has been measured (Roberts et al 



2004] ITrentham et al. 2005 Roberts et al. 2007; Blanton 



) - galaxies: dwarf - galaxies: luminosity function - 
fundamental parameters (classification, luminosities. 



with models if it is assumed the observed dwarf galaxies 
reside in relatively massive d ark matter halos wh ile lower 
mass halos remain invisible (Stoehr et al. 20021. If star 
formation efficiency is lower in lower density halos, or if 
large gas reservoirs are expelled from small potential wells 
after an initial burst of star formation and subsequent su- 
pernova activity, low mass galaxies will remain fairly dark 
and difficult to detect. Greater numbers of visible dwarfs 
might then be expected to inhabit the high density, high 
gas pressure environment of rich clusters if the intraclus- 
ter medium served to contain the dwarf galaxy outflows 
( Babul fc Rees|1992 ). For very small masses, cooling times 
may be greater than the age of the universe (Haiman et al. 
1996j ). Destructive processes which act in rich environ- 



et al. [2005 Barkhouse et al.| | ^ 

Several explanations to resolve this discrepancy between 
theory and observations have been proposed. Cosmolog- 
ical solutions include modifying t he power spectrum on 
small scales ( Kamionkowski fc L iddle 2000), perhaps by 
invoking warm dark matter TToh'n et al. 2000). Numeri- 
cal simulations generally predict subhalo counts as a func- 
tion of total mass or maximum circular velocity. For Lo- 
cal Group satellites, these values are not well determined, 
making direct comparison between these models and ob- 
servations difficult unless other paramete rs, such as the 
mass within a specific radius, are used (Strigari et al. 
20071. However, these more appropriate comparisons still 



p 

ments such as clusters may decrease dwarf number counts 
or lower luminosities making them more difficult to de- 
tect. Such processes includ e harassment, strangulation, 
and ram pressure stripping (M oore et al . 1996': 'Shay a fc] 
'Tully ||1984[|Bekki et al. 2001b,, Mori fc Burkert,20 00 ) aM 
would be expected to play a larger role in shaping galaxy 
populations in denser, more massive, and hotter clusters. 
Instead, high density regions have been observed to con- 
tain a larger fraction of dwarf to giant galaxies than low 
density regions. 

Another theory argues that reionization suppresses the 
collapse of gas in low mass halos thereby producin g a shal- 
lower faint-end slope in the luminosity function (Thoul fc 



find predictions inconsistent with observations. 

A large number of astrophysical solutions have been put 
forth to explain the apparent underabundance of dwarf 
galaxies. Number counts can be brought into agreement 



Weinberg 
tion in ha 



19961). This process will not affect gas accre- 
os that collapsed before the era of reionization 



- those that resided in denser regions, such as rich clus- 
ters, would preferentially collapse earlier when intergalac- 
tic gas was still cold and halos could accumulate the gas. 
Reionization will therefore predominantly suppress baryon 
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accumulation in low mass halos in lower density regions 
such as the poorer group and field environment, naturally 
explaining the trend of lo wer dw arf-to-giant ratios found 
in lower densit y regions dTully et al. 2002|). In support 
of this theory, Strigari et al.l (|2007p and ISimon & Geha 



(20071 find good agreement with theoretical expectations 



if the observed Local Group satellites correspond to the 
earliest collapsed halos, or the most massive at the time 
of accretion. 

Distinguishing between all these possibilities has proven 
difficult, but before we can draw any conclusions as to 
the source of the discrepancy, we must first constrain the 
number counts of dwarf galaxies and the faint-end slope 
of the luminosity function. No complete census of dwarf 
galaxies has been obtained in any environment, although 
much recent progress has been made in the detection of 
faint Local Group satellites. With SDSS survey releases, 
the number of known Milky Way and M31 satellites has 



nearl y doubled in just the past 4 years (Simon & Geha 
20071. This new information has mitigated the missing 



galaxy problem, but has not brought the number counts 
into alignment with ACDM predictions. Due to the sky 
coverage required for a complete survey and restrictions 
caused by the zone of obscuration, a complete census of 
Local Group dwarf galaxies is not on the horizon. 

The determination of the faint-end of the galaxy lumi- 
nosity function is hampered by the limitations in detecting 
the very faintest and lowest surface brightness galaxies, 
and this is compounded by the the great difficulty in de- 
termining accurate distances to, or ascertaining group and 
cluster membership for, such faint objects. To measure the 
faint-end of the luminosity function for nearby clusters of 
galaxies, researchers have often resorted to using statist! 
cal means to estimate cluster memberships (see e.g . |Bern- 



stein et al | ( |1995[ ); [TFentham| (|1997|); [Chiboucas fe Mat'eo 
( 2007p ). Others use colors or morpholo gy in an atternpt to 
directly assess memb ership probability ( Seeker et al.|1997 



Mahdavi et"ar]|2005[ ) . All these methods are fraught with 



large uncertainty. Spectroscopy can be used to establish 
membership, but to observe the large numbers of galaxies 
to the necessary depth is prohibitively time consuming. 

The M81 Group provides an alternative means to di- 
rectly establish membership for galaxies to magnitudes as 
faint as M/j ~ — 7 in a low density environment. Residing 
at a distance of only 3.6 Mpc, the galaxies in this group 
resolve into stars. Supergiants and giants near the tip 
of the red giant branch (TRGB) resolve over the diffuse 
main sequence population light. The degree of resolution 
is enough to establish the group membership status with 
high confidence. From color-magnitude diagrams, one can 
determine TRGB distances to unambiguously assign mem- 
bership. Given this fortuitous circumstance, we can probe 
the faint-end of the LF in a group environment down to 
Mr ^ —9 to high accuracy. 

The M81 Group is embedded in a filament of galaxies 
that runs up from behind the Galactic plane, shown in Fig- 
ure [T] The M81 Group has traditionally been considered 
to include the galaxies around NGC 2403 and those around 
NGC 4236, each about 800 kpc removed in projection from 
M81. Maffei 1-2 and IC 342 are considered to be the dom- 
inant members of a separate group, 2 Mpc away in projec- 
tion. The nearest neighbors at higher super galactic lon- 



gitude are a group around MlOl at SGL=64, SGB=-H23 
and the near Canes Venatici group at SGL=76, S GB=-|-9. 
It has been shown by Karachentsev et al. ( 2002a I that the 
subcomponents of the M81 Group around NGC 2403 and 
M81 itself are falling toward each other, consistent with a 
total mass within the zero velocity surface bounding the 
group of 1.6 X 1O^^M0. The interest of the present study is 
the more dynamically evolved subcomponent around M81 
contained within the box in Fig. [l] 

We surveyed a 8 x 8 degree area centered on the M81 
Group at a = 10:05:00 5 = 68:15:00 J2000 to search for 
new dwarf members of the group. Apart from the five 
brightest systems: M81, M82, NGC 2976, NGC 3077, 
and IC 2574, the inspected area contains 13 previously 
known dwarf galaxies which were discovered in the 1960- 
80s on POSS-I prints and the Tautenburg Schmidt tele- 
scope plates: Holm I, BK3N, KDG 61, Holm IX, Anon 
0952 (Arp's loop). Garland, BK5N, KDG 63, KDG 64, HS 
117, DDO 78, DDO 82, and BK6N 
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11985"). Later, 4 low surface brightness dwarfs: F8Dr, YMl, 
KK 77, and IK N were found on POS S-II films by, Caldwell 
et al. (1998).F roebrich fc Meusinger| (2000), Karachentsev 
et al.| ( |2000| ), and [Karachentsev et"al.| ( |2004| ) 



Here we report on the discovery of 22 new M81 Group 
galaxy candidates. In Section [2| we present the observa- 
tions and data reduction for our survey along with a dis- 
cussion of the dwarf galaxy detection techniques, survey 
limits, and follow-up observations. We provide measure- 
ments of the properties of the newly discovered dwarf can- 
didates in Section [3.1[ and construct the total group LF in 
Section 3.2 We discuss the dwarf population character- 
istics in Section 14.11 the distribution of M81 satellites in 



Section |4.2[ and implications of our survey results on the 
faint-end of the galaxy luminosity function and for cos- 
mological models in Sections |4.3[ and [4.4[ Conclusions are 
presented in Section [5j Throughout this work we assume 
a distance modulus to M81 of 27.80 (3.6 Mpc). 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Initial survey 

The M81 Group was chosen in large part because it is 
near enough at 3.6 Mpc that stars near the tip of the red 
giant branch just resolve in ground-based images, enabling 
derivation of the group luminosity function complete down 
to very faint magnitudes. Stars in Local Group galaxies 
resolve down to fainter magnitudes but a complete sur- 
vey would require all-sky coverage. For a comparable sur- 
vey of dwarf satellites around M31 alone, imaging of over 
800 square degrees would be required, including observa- 
tions within the zone of obscuration. The scale on which 
dwarfs cluster around parent galaxies is found to depend 
on the mass of the par ent halo accordin g to the relation 
r2t = 193(Mi2)^/^ kpc ( [Trentham et al.|[2006j where M12 
is the mass of the halo in units of IO^^A/q and is the sec- 
ond turnaround radius of spherical collapse, similar to the 
overdensity radius r2oo- At a distance of 3.6 M pc and halo 



mass of 1.6 x IO^^Mq ( Karachentsev et al.|[2002a| , M81 



has a second turnaround radius of 230 kpc, or only 3.6 de- 
grees. A survey of this group with complete areal coverage 
is therefore made possible with large format CCD detec- 
tors such as the square degree MegaCam detector mounted 
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on the 3.6m Canada-France-Hawaii Telescope. In addi- 
tion, the M81 Group is interesting in its own right as a 
rich group with a large late-type dwarf galaxy population, 
unlike our own well studied dSph-rich Local Group. 

MegaCam is a wide field camera composed of a mosaic 
of 36 CCDs that cover a total area of 0.96 x 0.94 deg with 
chip gaps 13 and 80 arcsec wide. Each CCD is composed 
of 2048 X 4612 pixels of 0.187". We imaged a 65 square de- 
gree region centered around M81 in a raster pattern that 
resulted in each location being observed twice for a to- 
tal integration time of 1096s. The survey was conducted 
entirely in the r'-band because we were primarily inter- 
ested in the resolved giant population. Typical seeing was 
0.75 arcsec, but ranged between 0.6 and 0.9 arcsec (Figure 

t. Most of the data were observed on photometric nights, 
though a few fields were imaged through thin cirrus. Ob- 
servations were carried out by CFHT staff in queue mode 
over 3 semesters in the period 2005B-2006. 

2.2. MegaCam data processing methods 

The raw images were initially processed by TERAPIX, 
an astronomical data reduction center located at the Insti- 
tut d'Astrophysique in Paris, which provides pipeline data 
reduction for MegaCam. The TERAPIX pipeline first 
evaluates data for image quality, identifies chip defects and 
cosmic rays, and generates weight and flag-map images. 
Each image is warped onto a common astronomical projec- 
tion (TAN) and coordinate system (J2000) using SWarp 
software and an astrometric reference catalog. SExtrac- 
tor is run on all images to perform the object detection 
and photometry. The astrometric and photometric cali- 
brations are determined using SCAMP. Astrometry is per- 
formed by downloading refere nce astrometic catalogs e.g. 
USNO-Bl dMonet et aL||2003[ ). These are cross correlated 
with each exposure and a pattern matching algorithm is 
run. Final astrometric uncertainties are dominated by er- 
rors in the reference catalogs. On large scales, the astro- 
metric accuracy is expected to be better than '--^ 0.3 arc- 
sec. Photometric calibration uses standard stars observed 
on photometric nights. The weighted, quadratic sum of 
SExtractor parameter MAG^UTO differences from over- 
lapping detections is used to put all fields onto a common 
zero-point scale. Residuals for bright objects are typically 
less than 0.05 magnitudes rms with even smaller zero-point 
uncertainties. These raw and processed MegaCam survey 
data are available online through the Canadian Astronomy 
Data Centre. 

Our fields were mosaicked in such a way that every re- 
gion was imaged twice with some additional overlap along 
the sides of each field. These pairs of 548 sec images had 
still to be stacked. Because MegaCam images are huge, we 
first split each image and its associated mask image into 
quadrants. The mask image contained fiagged chip defects 
such as bad columns and pixels as well as cosmic rays. 
Images were combined using IRAF tasks for image regis- 
tration and image combination. Pixels that were fiagged 
in each image were ignored during the average operation. 
Two to four individual images were used to produce final 
stacked quadrant images of size 30 x 30 arcmin. 

SExtractor was subsequently re-run on these combined 
images. Weight images were created during the stacking 
step with the weight of each pixel set to the square of the 



number of image pixels averaged together for the final pixel 
value. This weight image was used during the SExtractor 
detection of objects to put a significantly lower weight on 
those regions such as chip gaps where data came from only 
a single image as opposed to at least two images as is the 
case over the majority of the field. This prevented SEx- 
tractor from detecting excess noise in these regions while 
still reaching to ^ 0.75 magnitudes brighter than full ex- 
posure limits. SExtractor parameters were tuned to go 
deep enough to recover most of the mottles in the known 
resolved M81 Group galaxies by setting detection and de- 
blending levels to quite low values. Object detection was 
set to Icr with a minimum area of 3 contiguous pixels while 
the object deblending parameter NTHRESH was set to Icr 
and DEBLENDJVIINCONT to 5E-7. These low detection 
settings have the side effect of inducing SExtractor to re- 
cover excess noise, primarily in the halos of bright stars 
and galaxies, and occasionally in the numerous galactic 
cirrus patches obscuring this part of the sky. We err on 
the side of picking up spurious candidates rather than los- 
ing real candidates. 

We checked the reliability of the photometry by com- 
parison with stars of known magnitudes. |Perelmuter fc| 
Racine ( 1995 1 provide a catalog of i?, V, R magnitudes for 
all objects in the vicinity of M81. We transform the R 
magnitudes to the SDSS r' system using the rough tr ans- 
forma t ion r' = V - 0.84(V-R) -I- 0.13 provided in Fuk ugita| 
et al. (1996a). The average difference in r' (Perelmuter 

this work) for stars with 17 < r' < 19.5 is 0.01 with 
(7 ~ 0.11 (Figure [3|. Since this is only an approximate 
color transformation, we take this as strong support for 
the reliability of the photometry. 

2.3. Detection of 22 new candidate M81 companions 
The tip of the red giant branch lies at Mj 



-4.05, 

Up ^ 0.8 and using 
color transformations fo r metal poor population II stars 
(|2006|), we find Mr' tip 



with {V - I) 



tip 



1.5. For {V - R) 



from Jordi et al. 



-3.1 ± 0.1. 

With our 18 minute combined exposures, we reach limit- 
ing magnitude s for re solved stars of r' ~ 25.0 — 25.5. From 
Schlegel et al. (19981 maps, galactic reddening Ar' in this 
region is typically 0.1-0.2 but as high as 0.7 magnitudes. 
Because the M81 Group is so nearby at a distance of only 
3.6 Mpc ((m — M)o — 27.8), we can resolve stars down 
to limiting magnitudes of Mr' ^ —2.9 ± 0.5. The radial 
distance within the group will also affect the limiting ab- 
solute magnitude by ±0.15. Thus, stars at the tip of the 
red giant branch are just at the edge of detection in our 
MegaCam images while young main sequence and red su- 
pergiant stars, and intermediate age AGB stars are well 
resolved. 

Two different methods were used to search the Mega- 
Cam images for new M81 Group dwarf galaxy candidates. 
The first method used visual inspection of the images. Fi- 
nal images were examined by eye by two of the authors (IK 
and KC) for new candidates by identifying objects that 
had resolved and semi-resolved faint stars. Good candi- 
dates also had an associated diffuse low surface brightness 
(LSB) component, although this was not observed in every 
case. 17 new candidate M81 dwarf members were detected 
in this way. 

The second method used to discover new candidate 
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dwarfs implemented a 2-point correlation auto-detection 
routine. The auto-detection code calculates the 2-point 
correlation using the natural estimator for objects having 
magnitudes 22 < r' < 27 on scales of 1.1-6.4 arcsec (0.02 
- 0.11 kpc at the distance of M81) and within boxes of 
size ~ 40 arcsec (0.7 kpc). The code is set up to use any 
box or scale size but we were focusing on smaller galaxies 
with a concentration of resolved stars. Presumably large 
diffuse fields of stars could be missed, although the code 
did pick up a number of these as well. As the box marches 
across the image, the correlation values within the box are 
summed and any region with an excess value is flagged. 
This value was empirically chosen to be 80 based on re- 
covery of previously known and newly visually detected 
galaxies, and based on tests of the recovery of artificial 
galaxies (see below). The value chosen maximized the ra- 
tio of the rate of recovery of artificially added galaxies to 
noise detections, but does then necessarily impose mini- 
mum magnitude and concentration limits on our survey. 
All regions that were flagged were checked by eye. Typi- 
cally 40 — 200 such regions were found in each square de- 
gree MegaCam image. The great majority turned out to 
be noise detections along still remaining bad columns, or 
an excess of noise spikes around halos of bright stars and 
galaxies or within cirrus. Several detections in each image 
proved to be distant galaxy clusters. The remainder of 
the detections include the M81 candidates and previously 
known group members. 

The code recovered all previously known dwarfs and 
most of the ones found by eye in this paper, the 3 ex- 
ceptions being d0957+70, d0944-H69, and dl016-f69. For 
these 3, too few resolved points were detected by SExtrac- 
tor. It did pick out 2 dwarfs, dl009-H70 and d0946-h68, 
which have more of a diffuse LSB component than a re- 
solved stellar population and are thus more likely to be 
galaxies at somewhat larger distances. We retain all 5 
questionable objects in our full candidate list but although 
we still believe dl016-l-69 to be a good candidate, we con- 
sider d0957-|-70 to be a possible artifact and the other 3 
to be potentially more distant galaxies. 

The code also found a number of candidates not picked 
out by eye. A couple of these appear to have both stellar 
and diffuse LSB components: dl006-|-67 is the best of these 
candidates with a resolved stellar component contained 
within a diffuse low surface brightness region. d0934-f70 
is hidden behind cirrus but appears to have both a faint 
resolved stellar component and excess diffuse emission vis- 
ible over the cirrus. dl048-|-70 appears to have a diffuse 
component but only a semi-resolved stellar component. If 
real, this may perhaps be a more distant dwarf. A couple 
more have an excess of faint stellar-like detections but in 
these cases the objects do not appear to be concentrated 
nor have an LSB component and may prove to be distant 
galaxy clusters or chance concentrations of stars. Coordi- 
nates and basic parameters for a total of 22 new candidate 
dwarf galaxy group members are listed in Table [l] 

We display thumbnails of 3 candidate M81 Group blue 
compact dwarfs (BCDs) in Figure |4] along with a mosaic 
of another 9 good dwarf member candidates in Figure [5| 
Each galaxy image is 1.5 arcmin on a side. In Figures[6]-^ 
we display candidates which may turn out to lie slightly in 
the background of the M81 Group or simply be artifacts 



in the data, foreground clumps of stars, or distant galaxy 
clusters. Figure [8 displays a map of our 65 square degree 
survey region and the distribution of the new candidate 
sources along with previously known group member galax- 
ies. We show a few examples of galaxy clusters detected 
by our 2-point correlation routine in Figure [9] 

we display histograms of the number 
of the SExtractor detections within a ra- 



in Figure 10 
counts, N(TOflX 
dius of 30 arcsec centered on each of the new M81 candi- 
dates. For the 3 BCDs and d0944-|-71, another bright can- 
didate, we construct histograms for detections out to a ra- 
dius of 45 arcsec. A rough estimate of the fore/background 
object counts is obtained by measuring the number counts 
in a same sized region 2.5 arcmin west of each galaxy. 
These are displayed as shaded histograms overlaying the 
candidate population histograms. 

2.4. Limits of our Detection Methods 

We tested our automated method of detecting dwarf 
candidates through the 2-point correlation of resolved stel- 
lar populations using false stars. We wish to determine to 
what limiting stellar densities and magnitudes our galaxy 
detection routines are sensitive. To simulate artificial 
dwarf galaxies, large numbers of stars are added to real 
images with the IRAF package ARTDATA and convolved 
with the seeing of the field. Stars were distributed accord- 
ing to an exponential spatial density profile with an ef- 
fective radius, re, ranging from 4-90 arcsec (0.07-1.6 kpc). 
The stellar luminosity function has the form of a powe r 
law with power ~ 0.3 for the RGB (Ma karov et~al]|2006 1 . 
This value was allowed to vary randomly between 0.25 - 
0.35 in the simulations. We add between 100 stars and 
120,000 stars to each galaxy realization ranging in mag- 
nitude from r' ~ 30.0 at the faint end to 24.7 - 25.3 at 
the bright end. We take the tip of the RGB magnitude as 
our bright limit (24.7) but include up to 0.6 magnitudes of 
extinction. We also include a number of simulations with 
a bright limit up to 22. This allows us to simulate galax- 
ies with a range of stellar populations, from those having 
entirely old stellar populations to those also containing a 
young, bright component. 

The STARLIST task was used to generate a file suit- 
able for input to MKOBJECTS. Simulated stars were first 
convolved with the image seeing and Poisson noise was 
included before objects were added to the real images. 
Sets of 40 artificial galaxies were added at a time to single 
MegaCam quadrant images and run through our detection 
procedure. A total of 1900 galaxies were generated in this 
manner. In Figure |11[ we display a few examples of our 
simulated galaxies. SExtractor was re-run on these images 
and the 2-point correlation automated method was then 
used to search for these artificial dwarf galaxies. 

Although we expect the detectability of a galaxy to de- 
pend only on whether the dwarf galaxy has a significant 
resolved stellar population, we find that the detection rate 
and limiting depth of our survey depends on a number of 
factors. The primary factors are the total number of re- 
solved stars and the concentration of these stars. These 
correspond to the typical measured quantities of total in- 
tegrated magnitude and surface brightness of the galaxy. 
Results of 1200 artificial galaxy realizations generated with 
typical 0.64 - 0.75 arcsec seeing and with extinction values 
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< 0.2 are displayed in Figure [T2j We plot the integrated 
total magnitude and half-light radius for each simulated 
galaxy. Open points are those objects which landed on 
single exposure chip gap regions of the MegaCam field, or 
fell on bright stars or galaxies and are less likely to be 
recovered. Large black diamonds represent real candidate 
dwarf galaxies discovered in this survey. The dashed line 
displays the magnitude at which recovery drops to 90% 
within 2-magnitude bins, while the solid line represents 
the 50% completeness limit. We find that we are able 
to detect compact galaxies with < 11 arcsec down to 
r' = 20.5, or about Mr' — —7.3. For galaxies with Rg 

< 30 arcsec, we are able to reach r' = 19.0 (-8.8) with 
at least 50% completeness. For galaxies larger than Rg 
= 30 arcsec (~ 0.5 kpc) we are at least 50% complete to 
r' = 18 (-9.8). Assuming real dwarf galaxy sizes of Re < 1 
kpc, we expect our survey is 90% complete down to about 
r' = 18.0, or about -9.8 at the distance of M81. 



In Figure 13 we display the survey detection limits in 
the mag — /i plane. Symbols are as before. Lines of 
constant effective radius are labeled. A curved line is 
drawn that best bounds the region within which simulated 
galaxies are recovered. This is the traditional curve that 
describes isophotal limits for exponential profile galaxies 
(Allen & Shu 1979). The completeness of this survey de- 
pends on other factors, but the recovery boundary is ad- 
equately described by isophotal limits. At a given magni- 
tude, as one goes to fainter surface brightnesses, recovery 
fails because of increasingly low concentrations of stars 
and a failure of the eye or 2-point correlation methods to 
distinguish these from the background. At a given sur- 
face brightness, as one goes to fainter magnitudes, objects 
become smaller with decreasing stellar population sizes. 
Therefore, while not necessarily isophotal limits, the de- 
tection limits imposed on this survey can be translated to 
limiting angular sizes and surface brightness. 

However, the detection efficiency will also vary by loca- 
tion in these fields. Because our survey just reaches the 
tip of the RGB, any galaxy with a pure old population is 
detected only at brighter magnitudes and surface bright- 
nesses in regions of high extinction and fields with worse 
seeing. We test this with 240 artificial galaxies generated 
with 0.5 magnitudes of extinction. See Figure 14 On 
average, the 50% limit in the recovery is 0.3 magnitudes 
brighter than for the sample added with 0—0.2 magnitudes 
of extinction. Although the number statistics in this sam- 
ple is small, this suggests that the completeness limit in 
magnitude depends directly on the level of extinction and 
that our completeness limits may vary by as much as 0.6 
magnitudes throughout our survey region. Similarly, the 
depth of the group may affect completeness limits. Vari- 
ations in the distance with respect to M81 of ±300 kpc 
affect apparent magnitudes by ±0.18 mag. This is a very 
small effect but implies that the dwarfs on the near side 
of the group with pure old populations are more likely to 
exhibit detectable resolved stellar populations than those 
on the far side. 

We test the effect of poorer seeing with 400 galaxies 
that were simulated with 0.9 arcsec seeing. Results are 
presented in Figure 15 The solid line denotes the location 
where completeness drops to 50% with 0.9 arcsec seeing, 
while the dashed line represents the 50% completeness for 



seeing no worse than 0.75 arcsec. We find evidence that 
seeing has an effect on the recovery of resolved dwarfs, and 
this is most significant for compact galaxies (Rg < 20 arc- 
sec) where the completeness shifts 0.8 magnitudes bright- 
ward. This is to be expected. In these cases worse seeing 
blends the light of these small objects to produce a smooth, 
less resolved low surface brightness object. With the small 
number statistics, it is not clear that at larger Rg the 0.2 
arcsec worse seeing will affect recovery by the measured 0.4 
magnitudes, but the results do indicate that seeing has an 
affect on the detection rates of these objects. Fortunately, 
only ^ 19 combined images, or 11% of our survey area, 
will be affected by such poor seeing. 

To a lesser extent, one might expect the bright cirrus 
patches present throughout much of this region to affect 
detection limits through increased sky noise. From artifi- 
cial galaxy tests, however, we do not find this to signifi- 
cantly affect our recovery rates. Regions of high galactic 
stellar density can also limit detectability for the more dif- 
fuse dwarfs. At galactic latitudes ranging between 39 — 44 
deg this is not too much of a problem, but a few fields ex- 
hibited non-uniformly distributed excess numbers of stars 
throughout the images, increasing the confusion level. 
This will primarily affect fainter, lower surface brightness 
dwarfs with fewer total numbers of resolved stars spread 
over a larger area. 

Finally, two areal coverage issues will affect recovery. 
Any small dwarf lying behind bright saturated stars or 
galaxies will likely be missed. From false galaxy tests, 
where the galaxies are randomly placed within the im- 
ages, we find less than 4% of recoverable galaxies remain 
undetected because they landed on bright objects. This 
primarily affects only the smallest galaxies which can hide 
entirely behind brighter foreground objects, although a few 
more extended simulations were obscured by bright stars 
or galaxies. Also, due to the survey strategy of observ- 
ing each field twice by offsetting between images by half 
a field, all chip gap regions in the MegaCam field are ob- 
served for a total of only 10 minutes. With half the expo- 
sure time as the rest of the field, we expect the magnitude 
completeness limits in these regions to be about 0.75 mag- 
nitudes brighter than the above described survey limits. 
Because the detected dwarf galaxies in M81 have angular 
sizes generally greater than 16 arcsec, the vertical chips 
gaps of 13 arcsec should not affect recovery of any but the 
tiniest galaxies if by misfortune they should fall in these 
regions. There are however two large horizontal chip gaps 
which have single exposure widths of 80 arcsec each across 
the length of each MegaCam field, and this gap size is es- 
sentially doubled in our mosaic of two images per field. 
Since all our newly discovered dwarf galaxies have angu- 
lar extents smaller than 160 arcsec, it is therefore entirely 
possible for dwarf galaxies to be lost in these regions. The 
total area of these single exposure regions, including the 
area at the top and bottom of our full mosaic, is 9.7 square 
degrees, or ~ 15% of the total survey area. 

Based on these results, we find that we are nearly 100% 
complete to r' = 17 and 90% to r' = 18(Mr' = -9.8). 
At fainter magnitudes, a large number of factors act to 
diminish our ability to detect group galaxies, although the 
primary factors are simply the limiting magnitude and sur- 
face brightness to which we are sensitive. Factors such as 
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worse seeing, higher extinction, areal coverage by brighter 
objects, and single exposure regions serve to push these 
limits to slightly brighter magnitudes. 

2.5. Follow-up Spectroscopy 

Follow-up spectroscopic observations were carried out 
with FOCAS on the 8m Subaru Telescope and SCORPIO 
on the Special Astrophysical Observatory 6m BTA tele- 
scope. Spectra were obtained for our three BCD candi- 
dates along with several globular clusters in a couple of 
previously known M81 dwarf galaxies. 

Observations with Subaru were made 22-24 Nov 2006 
using the R300 / mm grating in second order and VPH 600- 
450 grating in first order. A two arcsec slit was used with 
the VPH grating, and images were binned 4x2, providing 
a spectral resolution of lOA at 4000A. For the R300 grat- 
ing, a 1 arcsec slit and 4x1 binning was used to obtain a 
spectral resolution of lOA. Spectral coverage ranged from 
3800 - 5250A for the VPH and 3650 - 5950A for the R300 
grating. Exposure times varied from object to object. 

Longslit spectra were obtained with FOCAS of globu- 
lar clusters in the previously known faint dwarf galaxies 
F8D1 and IKN. The globular cluster i n F8D1 was orig- 
inally noted by Caldwell et al. (19981. For two bright 
candidate galaxies (both BCDs) dl028-H70 and d0958-f 66, 
we obtained spectra by simply placing the slit across the 
long axis of the galaxy. These data were processed us- 
ing standard IRAF routines. Images were first cosmic ray 
corrected using CRMEDIAN and then bias subtracted us- 
ing 10 median combined bias frames. Arcs were observed 
before or after each observation and used to wavelength 
calibrate the spectra. Lines in the arcs were identified 
using IDENTIFY and REIDENTIFY, and FITCOORD 
was used to obtain a dispersion correction solution which 
was applied using TRANSFORM. The two-dimensional 
spectrum for each object in the sli t was extracted using 
APALL. Finally, RVSAO/XCSAO ( [Kurtz fc Mink|[l998 l 
was used to obtain a cross correlation solution from tem- 
plate spectra for the redshift of each object and hence the 
galaxy itself. Spectra are displayed in Figure 16 Radial 
velocities are listed in Tabled 

Further observations were performed for the 3 BCD can- 
didates with SCORPIO in long-slit mode installed at the 
BTA 6-meter telescope prime focus. The grism VPHG400 
was used with a 1 arcsec slit yielding the spectral range 
~ 3500 — 7500 A chosen to include Ha with a resolution of 3 
A/pixel and FWHM - 20 A/pixel. The objects dl028-h70 
and d0958-t-66 were observed during the night 11 Novem- 
ber 2006 with the exposure time of 600 sec under 2" seeing. 
The galaxy UGC 5497 (dl012-h64) was observed on 20 Oc- 
tober 2007 with an exposure time of 900 sec under 2.5" 
seeing. Spectrophotometric standard stars were also ob- 
served for calibration. The standard pipeline with the us e 
of MIDAS, as described in [Afanasiev fc Moiseev] ( |2005| ), 
was applied for the reduction of these long-slit spectra. 
All three BCD candidates show knotty Ha emission. We 
display the spectra in Figures [iT] - [19] Due to the short 
exposure times and low S/N it is difficult to measure line 
ratios. However, we note that the [Nil] A6584 line is quite 
weak in comparison to Ha, especially for d0958+66. We 
estimate an upper limit for the line ratio of 0.05 fo r this 
galaxy. From the relations of Kewley fc Dopita ( 2002, ) and 



Pettini fc Pagel ( 2004 ) , this leads to an upper limit for the 
metallicity of 0.3 Zq. Mean heliocentric radial velocities 
of the objects are included in Table [2] 

2.6. Multi-band Imaging Follow-up 

Using the UH 2.2m with Tek 2K CCD on 19-24 Dec 
2006, we imaged 9 of the dwarf galaxy candidates in V and 
/—bands. Seeing ranged from 0.7— 1.4 arcsec and 5 nights 
were photometric. Images were 4 x 900s in the F— band 
and 5 x 720s in the /—band. Standard stars were observed 
on each photometric night. Images were processed using 
standard IRAF routines. Due to the slow speed of the 
mechanical shutter, a shutter mapping correction was also 
implemented. This was necessary for the short exposure 
images (5-20s) of the standard stars. 

In most cases, the images did not go deep enough to 
produce useful color-magnitude diagrams (CMD). For 8 
galaxies observed on photometric nights, we measure in- 
tegrated V — I colors within circular apertures out t o th e 
effective radius as measured in the r' band (Section [3. 1[). 
We co rrect for extinction using the maps of [Schlegel et al. 
( 1998 1. V — I colors are hsted in Table [s] We also gener- 



ate color images of the 9 candidates. To create images for 
the blue channel we take 2 *Vi 



image 



limage- The V and 



/-band images are taken to be green and red respectively. 
Color images are displayed in Figure [20] 

Large young, blue stellar populations are evident in 
d0926-l-70, d0958-f-66, d0959-l-68, and dl028-l-70. These 
are objects classified by us as BCDs, dis, and tidal dwarfs. 
These multi-band data may not be going deep enough to 
reach the tip of the red giant branch. While resolved stars 
are seen in the MegaCam survey images, the lack of visible 
resolved stars in the UH 2.2m images of d0955-|-70 suggests 
it may consist predominantly of an older stellar popula- 
tion. Similarly if members of the M81 Group, d0934-h70 
and dl006-|-67 also appear to contain predominantly older 
stellar populations. 

We have recently obtained HST/WFPC2 multiband 
imaging for 14 of the candidates and have been awarded 
time to complete observations of the sample with 
HST/ACS. These data will be deep enough to measure 
distances from the CMD tip of the red giant branch. Re- 
sults will be presented in a subsequent paper. 

3. RESULTS 

3.1. Properties of the M81 companion candidates 
To determine centroids of the new objects, we used the 



around 



intensity weighted means of pixel values, ^ j 
each galaxy. EUipticity and position angle around these 
centroids were th en derived from the in tensity weighted 
second moments (Seeker fc Harris 19971. The measured 
ellipticities, (1— 5/a), for the faintest objects are highly un- 
certain due to low S/N of the diffuse light component and 
small visible stellar populations. Measuring total magni- 
tudes for these very faint and diffuse extended resolved 
galaxies is challenging because of expected contamination 
from foreground stars, and in several cases, from excess 
foreground cirrus reflected flux, not to mention extinc- 
tion. We therefore estimated total magnitudes using sev- 
eral methods. 

We first mask all obvious bright foreground stars and 
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background galaxies around each of our candidates. Sim- 
ple aperture magnitudes were obtained for each galaxy 
with radii 50, 100, 150, or 200 pixels (9.35, 18.7, 28.05, and 
37.4 arcsec) depending on the visible extent of the galaxy. 
The sky was determined in an annulus 4 pixels wide with 
a radius of 400 pixels. Alternatively, we measure total flux 
in an identical sized aperture in 4-5 regions around each 
galaxy. Taking the average of these as an estimate for the 
contribution from both sky and foreground/background 
object flux we subtract this from our object flux for a sec- 
ond estimate of total magnitude. 

To obtain structural parameter measurements for each 
galaxy, we extracted fluxes in successive elliptical aper- 
tures around each galaxy out to a radius of 40 arcsec and 
truncated the data if it exhibited a > 30% increase in slope 
attributed to contamination by an unmasked neighboring 
object, or, more often, to features like bright cirrus bands. 
The central 4 arcsec were also discarded. We then flt both 
the surface brightness profile with a Sersic function and 
the curve of growth with a cumulative Sersic function. We 
extracted parameters from the best of these two fits. This 
provides a third estimate of the total magnitude along with 
central surface brightness, effective surface brightness, and 
half-light radius. The generalized Sersic function is defined 
as 



Hr) = lo e 



-(r/r„)i/" 



(1) 



where n is the Sersic parameter. For n = 1, this reduces 
to an exponential profile which usually provides a good fit 
for disk and dwarf galaxies, and for n = 4.0 becomes a 
de Vaucouleurs profile. The cumulative Sersic function is 
then 

I{r) = 2Traorln ■ j[2n, (r/ro)^/"] (2) 
where 7[a,x] is the Incomplete Gamma function, 

exp{-t)e-^dt. (3) 



This can be integrated to obtain the total object fiux. 

It = 2TTaorlnY[2n] (4) 
with total magnitude, 

mt = iio- 51og(ro) - 2.51og(nr[2n]) - 2.0 (5) 



(Binggeli & Jerjen 1998) where T[a] is the Gamma func- 
tion. 

The fitting of this nonlinear function to t he data was 
done using a Levenberg-Marquard algorithm ( Press et al. 
1992 1 which performs a minimization that implements 



an inverse-Hessian method far from the minimum and 
switches to a steepest decent method as the minimum is 
approached. In Figures [2T] - [26] we display both curve of 
growth and surface brightness profile fits for our candidate 
galaxies. From these fits we directly obtained the central 
surface brightness, scale length, and profile type, where 
scale length in the generalized Sersic function is directly 
related to the half-light or effective radius through 

i?y" = i?y'V(2.3026&„) (6) 
with b„ = 0.868 n - 0.142 ( |Caon et al.|[l993 l. 



The majority of the galaxies are best fit with Sersic in- 
dex n < 1.0. From inspection of the plots it is clear that 
these low n fits are primarily driven by large cores, typical 
of dwarf galaxies. In the outer radii, the surface bright- 
ness profile tends to drop exponentially. We overlay best 
fit exponential profiles to the outer radii in the surface 
brightness profiles of Figures [2l] - |26j A number of cases 
exhibit a sharper drop-off, but arc still consistent with an 
exponential fall-off within the large uncertainties of the 
low surface brightnesses at these radii. 

We estimate typical magnitude measurement errors of 
0.5. Sources of uncertainty come from sky determination 
(0.4, estimated from curve of growth fitting when using 3 
different estimates of sky values), foreground/background 
contamination (typically 0.2, from magnitudes measured 
with and without masking of likely foreground sources), 
aperture magnitude measurement error from qphot (0.04), 
and the zerop oint uncertainty (0.015 , estim ated from com- 
parison with |Perelmuter fc Racine (19951 catalog stars). 
Cirrus may also effect measurements of some of the galax- 
ies. Sky estimates in these regions have an even larger 
uncertainty, and faint old stellar populations may be lost 
in the excess sky noise. In particular, d0934-|-70, if real, 
appears largely obscured by cirrus. d0944-|-71, d0959-|-68, 
and d 1048-1-70 are strongly effected by fore or background 
objects. We estimate total magnitude uncertainties of 0.8 
for d0958-|-68 and 0.6 for dl048-h70. The magnitude for 
d0944-f71 is corrected by assuming a symmetric shape. 
Flux from the star superposed quadrant is replaced by the 
ffux in the symmetrically opposite quadrant. We assume 
a magnitude uncertainty of 0.8. Candidate dl019-f69 is 
affected by missing data in a chip gap; the measured mag- 
nitude is corrected in a similar way as d0944-f71. 

We estimate measurement errors in Rg by compar- 
ing results from the surface brightness profile and curve 
of growth fits and comparing differences from fits using 
slightly different estimates for the mean sky value. We 
find fractional uncertainties in Re range from 1% for the 
bright BCD dl028-|-70 to 40% for low surface brightness 
object d0955-|-70, with a median uncertainty of 10%. Sim- 
ilarly, we find an uncertainty in the central surface bright- 
ness ranging from 0.05 to 0.5 magnitudes arcsec"^ with a 
median uncertainty of 0.2. 

Because of the shallow nature of our survey, we can as- 
sume that our magnitude and size measurements are lower 
limits. To obtain estimates of total magnitudes, all mea- 
sured magnitudes had to be corrected for light lost in the 
sky noise. Since many of the profile fits were infiuenced by 
the bright 0/B stars, it is possible that measured profiles 
do not accurately refiect that of the underlying old stellar 
populations. We therefore calculate the missing fraction 
of light assuming both the fitted profile shape and an ex- 
ponential surface brightness profile. Dwarf galaxies may 
be better fit with King profiles than exponential due to 
the presence of large cores and possible truncation in the 
wings of the profile. This is supported by our generalized 
Sersic function profile fits which find n < 1.0 for most of 
our galaxies, implying a steeper fall off of the galaxy light 
in the wings of the profile. Magnitude corrections calcu- 
lated by assuming exponential profiles may therefore be 
overestimates for these galaxies, and thus provide upper 
limits on the total galaxy luminosities. 
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From an isophotal hmiting surface brightness (ranging 
from 27.0 to 28.5 depending on the sky level and presence 
of cirrus) , along with the measured central surface bright- 
ness and total magnitudes, we extrapolate to obtain the 
intrinsic total magnitude. Over x scale lengths, x — r/ro, 
the surface brightness of an exponential profile drops by 



2.5 log e" 



1.086a; 



for 71 = 1 



1.086a;i/" 



(7) 
(8) 



for a generalized Sersic profile. The number of scale 
lengths observed between /io and /iHm is 

^x^{^iu^~ ^io)/lm& (9) 



or 



Ax= ((/^i,™-Mo)/1.086)". 



(10) 



For niT - = rriextrap with 

TTii oc 2.5 log[(l + x)e~^] (exponential) (11) 

and 

mj(x2.51og[/ e""" "rdr] (Sersic), (12) 
Jo 

we find the extrapolation beyond the measured magni- 
tudes is 

Am = 2.5 log[l - (1 + Ax)e-'^'=] (13) 



and 



Am = 2.51og[7[2ri, Ax^/"]/r[2n]] 



(14) 



for exponential and generalized Sersic profiles respectively. 
See also TuUy ct al. (1996). Exponential profile corrections 
from the isophotal magnitude are less than —0.1 mag for 
our brightest new candidates, but as high as —1.2 mag for 
the faintest objects with a mean value of —0.4. Correc- 
tions for the best fit Sersic profiles range from 0. to —0.8 
magnitudes with a mean of —0.2. Besides assuming a pro- 
file shape, these corrections further assume an old stellar 
population exists which may be too faint to detect, even 
as a diffuse low surface brightness component, in some 
cases. For potential tidal dwarfs like d0959-|-68, it is not 
clear that this assumption is valid. There may be evidence 
that d0959-|-68 is dominated b y a very young popu lation 
formed less than 70 Myr ago (Durrell et al. 2004 1. We 
therefore do not apply a correction for this object. 

We provide the detailed structural and photometric 
properties of these new objects in Table [3] If we assume 
all candidates are group members at the distance of M81, 
we find Mr' magnitudes range from -6.7 for d0944-|-69 to 
-13.3 for dl012-|-64. Half-hght radii range from 90 to 490 
pc for objects well fit with Sersic functions. Because we 
do not yet have distances to these galaxies, we estimate 
that intrinsic magnitudes and radii can vary by as much 
as 0.14 mag and 6% in size (assuming a distance range 
of ± T2t)- We provide r' total magnitudes of previously 
known M81 Group members that lie within our survey 
region in Table |4] For the largest galaxies, we expect 
larger magnitude measurement errors, up to 1 mag, due to 
large errors in the local sky determination and subsequent 
negative features in the images. For measuring the local 
sky backgrounds, the TERAPIX reduction pipeline imple- 
ments box sizes much smaller than the size of our largest 
galaxies. 



3.2. M81 Group Luminosity Function 

To construct the M81 Group differential luminosity 
function (LF), we bin number counts as a function of mag- 
nitude in 2-magnitude bins. Such large bin sizes are nec- 
essary because of the small number counts we are working 
with and because of large magnitude measurement errors 
(> 0.5 mag). Counts are normalized by the size of our 
survey region. We construct luminosity functions for all 
previously known members, all previously known plus all 
new candidates from this work, and for previously known 
plus best candidates from this work. In the case of 'best', 
we include objects in Figures 4][5 We have excluded counts 
for objects we believe may be background or which could 
potentially prove to be artifacts. Of the previously known 
objects, we also exclude from this 'best' sample Arp's loop 
and Garland since we do not have reliable magnitude esti- 
mates for these likely tidal objects. The luminosity func- 
tions for these different samples are displayed in Figure |27j 
We have constructed the LFs using measured magnitudes 
(top) , exponential profile corrected total magnitudes (mid- 
dle), and Sersic profile corrected total magnitudes (bot- 
tom). 

Based on our artificial galaxy tests, we shade the re- 

17, 

18, to be 



gion where incompleteness sets in. Brighter than 
we expect we are 100% complete, and to 



at least 90% complete. We may still be missing a small 
number of dwarfs which are lurking behind bright stars or 
other galaxies, but find this should be less than a 4% ef- 
fect, primarily affecting only the smallest objects. We use 
our simulations to determine estimates for our complete- 
ness at each magnitude bin in the LF. The lower limits in 
the detection efficiency imply that, at least in regions with 
good seeing and low extinction, we should recover all ar- 
tificial galaxies down to r' — 17. At magnitudes brighter 
than this, we find no candidate dwarf galaxies with ef- 
fective radii greater than 30 arcsec. Fainter galaxies will 
presumably be smaller. As a conservative estimate, we 
assume a maximum Rg of 35 arcsec and uniform distribu- 
tion in Re and calculate the fraction of galaxies we would 
expect to recover in each 2-magnitude bin of our luminos- 
ity function. In reality, we expect galaxies to be weighted 
towards smaller sizes at fainter magnitudes. The assump- 
tion of a uniform distribution should therefore provide us 
with an upper limit for the fraction of the real population 
we would expect to have missed, and an upper limit for 
corrected counts. Results are plotted in Figure [28] The 
solid line is a fit to the recovery results of simulations with 
good seeing, low extinction, and falling on full integra- 
tion regions. The dashed line incorporates the 15% of the 
survey with single exposures and takes into account any 
lower completeness due to worse seeing and extinction in 
each field. We apply these fractions as completeness cor- 
rections to the counts in our luminosity function (open 
circles in Figure 27 1. Applied corrections also include a 
3% correction for areal coverage by brighter objects for 
r' > 17.0. Although these are only crude corrections since 
we do not know the true distribution of M81 Group galaxy 
sizes faintward of our 100% completeness limits, we do not 
find any evidence suggesting the faint-end slope steepens 
beyond what we measure from our observed galaxy counts. 

The lumino sity distribution of galaxies is traditionally 
modeled by a Schechter ( 1976 ) function or with a simple 
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power law for the faint end only. To measure the faint- 
end slope of the magnitude distribution, we perform a 
power law fit to the complete set of galaxy counts be- 
tween 8 < r' < 18 and find a faint-end slope, a, ranging 
from — 1.26j;QQ8 f™' the case of measured magnitudes to a 



maximum of — 1.39_qq5 for exponential profile corrected 
total magnitudes. The slope for the full, completeness cor- 
rected sample ranges from -1.28j;!^:^5 -l-32to.05- We 
measure a slightly shallower faint-end slope for only previ- 



ously known group members of — 1.191q to 



4.26: 



-0.08 
0.08- 



However, it is clear from the figures that the slope does not 
increase much with the addition of these new candidate 
dwarf group members; only the completeness has shifted 
to fainter limiting magnitudes. Completeness corrected 
counts as a function of both measured and Sersic corrected 
magnitudes are well described by a — —1.3 slopes. It is 
possible that our luminosity functions display an upturn at 
the faint-end, at magnitudes r' > —15. We therefore also 
fit all counts at the faint-end in the range —15 < r' < —10 
where we expect to be nearly complete. This provides us 
with the steepest slopes consistent with our data, with 
a ~ —1.45. However, we caution that this measurement 
is based on few counts and low number statistics. Best fit 
values for the faint-end slope, a, for all samples are listed 
in Table E 

Due to the overall small number counts of group galax- 
ies, uncertainties in each bin of our differential luminosity 
function are large. This has the effect of producing a num- 
ber count distribution that does not monotonically rise as 
a function of magnitude and, with so few data points, is 
not well fit with a power law. We therefore also produce 
a plot of the unbinned cumulative distribution of galaxy 
counts as a function of measured and corrected total mag- 
nitude for all subsamples in Figure 29 We fit a cumulative 
Schechter function to the cumulative counts. 



$(M)rfM = 0.41nlO(^, / ioO-4("+i)(M.-Af)g-io' 



r)0.4(M* -M) 



such that 



(15) 



iV(< M) = cl3,-f[a + 1, iqO-MM.-m)^ 



(16) 

using maximum likelihood techniques with a Poisson es- 
timator, is not only not well constrained due to the 
very small number of objects at the bright end of the LF, 
but is also affected by large errors in measured magnitude 
for the the brightest galaxies. For M81, we know our mea- 
sured total magnitude is strongly affected by improper sky 
subtraction. For the luminosity function we therefore as- 
sume a total magnitude for this galaxy, M^' = —21.8, from 



functions of 'best' and 'all' candidate samples with Sersic 
corrected magnitudes. It can be seen that M* is uncon- 
strained at the bright end, but the choice of Af* will not 
strongly effect the value of a. 

Overall, we obtain very similar values for the faint-end 
slope as compared with the power law fit to the faint end 
of the differential LF. For all samples, we find a faint-end 
slope a ~ —1.28. Slopes for only originally known and 
'best' candidates are insignificantly shallower than for all 
candidates. We furthermore find the faint-end slope to be 
robust against corrections from measured to total magni- 
tude. 

Fits to completeness corrected counts do not produce 
steeper faint-end slopes, but do extend our fits to all counts 
down to r' ^ —7. Since we expect these completeness 
corrections are already an upper limit given the expected 
true size distribution of real galaxies, we do not believe 
the slope steepens appreciably at magnitudes fainter than 
our completeness limit. The steepest measured slopes are 
for the full sample of objects which assume all candidate 
objects we have detected are real galaxies within the M81 
Group. Certainly some subset of these candidates are not 
true M81 Group galaxies, so we can take these slopes as an 
upper limit. We find the best fit slope for the subsample 
of most likely M81 members using Sersic profile corrected 
magnitudes to be -1.30 ± 0.10 and -1.26 ± 0.05 for the 
differential and integrated LFs respectively. 

4. DISCUSSION 

4.1. Dwarf Population Characteristics 

Of the new candidate galaxies, roughly 45% are believed 
to belong to the dwarf spheroidal class with dominant old 
stellar populations, while the rest exhibit dwarf irregular 
or BCD morphology with younger stellar populations and 
ongoing/recent star formation. At least one new object 
^^^may be a tidal dwarf, formed in the tidal streams connect- 
ing the three closely interacting galaxies of M81, M82, and 
NGC 3077. This area contains a number of knots of star 
formation and tidal debris. Previously known or suspected 
tidal objects include Arp's loop knots. Garland knots. Ho 



IX, and BK3N ( [Makarova et al.„2002| . These objects are 
found to be dominated by young stellar populations and 
have high gas contents (de Mello et al. 2008 Sabbi et al. 
|2008 ). Davidgc (2008| recently found an arc of stars in the 
M82 halo which may have formed around the same time 
as these other tidal features. One of our candidate ob- 



previously obtained R band data (Tully & Pierce 20001. 
For fitting the cumulative Schechter function to our data, 
we hold M* and <?!>* constant while solving for the best fit 
slope, a. Results for the different samples are presented 
in Table [6) To test the effect of on the fit, we re- 
fit the data allowing all parameters to vary. We find the 
slope changes typically by about 0.015, and no more than 
0.04, smaller than the random errors in our fitted slope 
values. Results are similar when we use our original total 
magnitude measurement for M81. Therefore, the value of 
Af*, although coupled to a, does not significantly affect 
the determination of the faint-end slope. In Figure |30] 
we display the la a — error ellipses for the luminosity 



jects, d0959-H68, was first noted by |Durrellet al | ( |2004[ ) in 

a search for intragroup stars in the field of the M81-NGC 
3077 HI tidal tail. It has recently undergone a burst of star 
formation as evidenced b y th e large population o f bright 
blue stars seen in Figure 20 Durrell et al. (20041 match 
isochrones to the CMD oFthis object and find evidence 
that star formation occurred 30-70 Myr ago, after the for- 
mation of the HI tidal arm. In a search for HI clouds 



in the M81 Group, Brinks et al. (20071 find a number 
of free floating clouds without optical counterparts which 
they suggest may be tidal debris from the 3-body inter- 
action involving M81, M82, and NGC 3007. One of these 
clouds is just offset from d0959-|-68, about 1'.5 arcmin to 
the SW, several times the measured optical half-light ra- 
dius. This HI cloudlet has a radial velocity distinct from 
both the Milky Way and M81 tidal arm. If it turns out the 
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optical component has a similar radial velocity as the HI 
cloud, it would be likely that the two were associated. It 
may be possible tha t, like the Local Group Phoenix dwarf 
(Young et al. 20071, blow out from supernovae energetics 
trom this young stellar population has displaced the neu- 
tral hydrogen gas. We display this new object along with 
previously sus pected tidal dwarf galaxies superimposed on 
the HI map of Yun et aL|pJ9 4) in Figure [31] 



Three new M81 candidate member galaxies (dl028+70, 
d0958+66, dl012+64) are of the poorly understood blue 
compact dwarf class. Spectra of these galaxies were ob- 
tained with Subaru/FOCAS and BTA/SCORPIO. Radial 
velocities of ~ -1-60, —100, and -1-150 km s~^ support group 
membership for d0958-f66, dl028-h70 and dl012-|-64 re- 
spectively. All spectra exhibit strong Balmer absorption 
lines indicative of recent star formation. The BTA 6m 
spectra also find Ha emission in all three cases indicat- 
ing active star formation, in agreement with the blue 
(F — /) ~ 0.8 colors found for these galaxies. These three 
candidates bring the total number of M81 Group BCDs to 
at least four, including UGC 6456 , long considered an iso- 
lated mem ber of the M81 Group (jLynds et ar||1998| [Tull yl 
et al.|198T l. DDO 82 may have a BCD component as well. 



It has an optical structure and color that corresponds to 
that of BCDs but has an HI flux contaminated by local 
HI. UGC 6456 lies far outside our survey region, at SGL 
= 37°, SGB = -hll° in Figure [T} while DDO 82 hes at the 
projected distance of the second turnaround radius for the 
M81 Group. BCDs tend to be found more frequently in 
field environments than clusters and are one of the least 



clustered galaxy types known (van Zee et al. 20011. The 
three BCDs discovered here all lie in low density regions at 
a minimum projected distance of 140 kpc from M81, with 
dl012-|-64 lying outside of the second turnaround radius 
at a projected distance of 315 kpc. The excess of this type 
of galaxy observed in this group may indicate that the ter- 
mination of star formation is due to processes such as ram 
pressure stripping which occur preferentially in more mas- 
sive and more dynamically evolved groups/clusters. How- 
ever, given the projected distances of the M81 BCD pop- 
ulation, it is possible that these galaxies are on first infall 
orbits and processes such as strangulation that serve to 



shut down star formation even in poor groups (Kawata & 
|Mu lchaey"2008) have yet to act on these galaxies. 



detection limit would enclose only the brighter and larger 
UCDs. Due to the extremely compact nature of UCDs, 
it is probable that the smaller ones would not appear in 
our images as resolved concentrations of stars. Rather, the 
hght on these small scales would be blended making de- 
tection based on a resolved stellar population difficult and 
in some cases impossible. However, we would expect to be 
able to detect the larger and brighter subpopulation. As 
we do not detect a single UCD, this suggests that UCDs 
may be absent in this poor, unevolved environment. 

No new bright (r' < —12) ultra-diffuse dwarfs like F8D1 
and IKN have been found in this survey. This could par- 
tially be due to selection effects since neither our eye nor 
the two-point correlation code would be particularly sen- 
sitive to extremely diffuse objects, although these two ob- 
jects were easily picked out both by eye and recovered with 
the two-point correlation technique. We would therefore 
expect to detect similarly low density objects if they ex- 
isted. We do find a number of very diffuse, faint, and small 
candidates including dl013-F68, dl014-h68, and dl016-h69 
which may or may not prove to be real objects. 

In comparison with the Local Group, we find a very 
different morphological mix of dwarf satellites. In the 
M81 Group, we have found 3 new BCD candidates while 
no such objects have been discovered around M31 or the 
Milky Way. The fraction of late type, dwarf irregulars is 
also much higher in the M81 Group. Including the recent 
detections of ultra-faint dwarfs in the Local Group, the 
fraction of late types is ~ 25% while in the M81 Group, in- 
cluding our new best candidates, the fraction of late types 
could be as high as 55%. This different morphological 
mix had been noted earlier for giant galaxies and brighter 
dwarfs, but recently discovered faint dwarfs in both groups 
maintain the distinctly different population fractions. It is 
inferred that the M81 Group is dynamically less evolved 
than the Local Group. 

We do not find any objects as faint as the recently dis- 
covered ultra-faint Local Group galaxies. This is due to 
incompleteness at these low magnitudes. We do find one 
candidate, d0944-|-69, with an absolute r' magnitude of 
—6.7 ± 0.5. This would make it the faintest known M81 
Group dwarf. It also would have one of the smallest sizes 
for known dwarfs of Re = 90 pc, similar in size to some of 
the recent Local Group satellite discoveries with Re rang- 



We find no ultra-compact dwarfs (UCDs) such as those ing from 23 - 125 p c (,Willman et al. 2005 



found in the riche r Fornax and Virgo clusters ( Drinkwater 
et al ||2000[J2003[ ). Neither have UCDs been discovered in 



the Local Group. This may support the hypothesis that 
UCD formation is driven by a mechanism that occurs pref- 
erentially in denser environm ents, such as tidal s tripping 
of nucleated dwarf ellipticals (Bekki et al.||2001al or late- 
type spirals ( Moore et al.|1998 1. Alternatively, Tellhauer fe| 



2006: Ma rtin et al.||2006[ [Belokurov et al.||2007a| |WalsE 



Zucker et al. 



et al.||2007 ). We do not yet know if this object is either a 
real galaxy or a member of the M81 Group. If this object 
does turn out to be a group member, its lower limit half- 
Hght radius would impinge on the 40 < Re < 100 pc 'size 
gap' b etween dwarfs and globular clusters ( Belokurov et al. 



Kroupa ( 2002 ) suggest the possible formation of UCDs as 
objects evolved from young massive star clusters which are 
formed within tidal tails of massive galaxy-galaxy merg- 
ers. In this scenario, one might expect the formation of 
such objects to be ongoing in the tidal tails of the galaxy 
interactions occurring in the M81 Group although we find 
no evidence for this. It is possible that we would be un- 
able to detect the more compact UCDs in this survey using 



our current search techniques. We show in Figure 13 the 
location in r' — fig space the UCDs would occupy. Our 



2007a I . However, as most of the galaxy light likely fails 
below the level of the sky, we are certainly underestimat- 
ing its physical size and perhaps, even with the magnitude 
correction, total luminosity. 

The other candidates all have properties consistent with 
normal faint dwarfs. We plot M^' vs log (half-light radius) 
in Figure [32j New candidates are shown as filled circles 
(assuming a distance modulus of 27.8) while previously 
known group members are displayed as filled squares. The 
three candidate BCDs are distinguished by their high sur- 
face brightnesses. For comparison, we include Local Group 
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dwarfs in this plot. Recent discoveries are denoted by open 
circles while previously known galaxies are shown as open 
squares. To transform fro m V to r' magnitudes, we use 
r' = l/-0.84(l/-i?)+0.13 ( jFukugita et al.|1996b[ ) assum- 
ing an approximate {V — R) color of 0.6 for dwarf galaxies. 
Lines of constant surface brightness are shown. New Lo- 
cal Group detections only became possible with the SDSS 
survey reaching effective surface brightne sses below 27.0 



magy arcsec (Belokurov et al. 2007b I, corresponding 



to a stellar surface density limit. Incompleteness in our 
survey may be starting to set in by 26.0 mag^' arcsec^^ 
and most of our detections are brighter than 27.0 mag^' 
arcsec"^, so we may still be missing a large population 
of ultra-faint dwarfs. According to Figures 13 and 32 we 



would expect these missing galaxies to have total magni- 
tudes fainter than Mr' ~ —8. 

4.2. Distribution of M81 Galaxies 
The asymmetric distribution of dSphs noted by 



Karachentsev et al. (2000j) has become less pronounced 
with the addition of the new objects. As discussed in that 
paper, a possible cause for the asymmetry is contamina- 
tion of the M81 Group area with many reflecting nebulae 
(galactic cirrus) which made it difficult to find extreme 
LSB dwarfs that are unresolved into stars or have only old 
stellar populations. In Figure [33| we show the M81 region 



of the sky in the Schlegel et al. (11998!) galactic dust map 



Superimposed are the locations ofgalaxies in the group. A 
number of new candidates (diamonds) can be seen to lie in 
regions of higher extinction and, given the dearth of detec- 
tions in regions with the highest extinction, it is possible 
we are still missing faint M81 dwarfs concealed by these 
clouds. However, it seems unlikely we would be missing 
dwarfs of comparable brightness to the ones displaying the 
noted asymmetry, given that the new candidates from this 
survey which are detected in regions of high extinction are 
all 1-2 mags fainter. 

We compare the projected distribution of morphological 
types in Figure |34j Late types, early types, and tidal can- 
didates are distinguished by symbol types while point sizes 
indicate surface brightness with brighter surface bright- 
nesses denoted by larger points. All the bright giant galax- 
ies are late types (circles) and all candidate tidal dwarfs 
(squares) lie near the center of the cluster right around 
M81, M82, and NGC 3077. Early type dwarf spheroidals 
(hexagons) lie primarily towards the group core while all 
late type dwarfs lie towards the periphery of the group. 
A clear morphology-density relation is seen. We find the 
median distance from M81 for late type dwarfs is 1.7 times 
that for dwarf spheroidals. 

4.3. Faint-end of the Group Environment Luminosity 

Function 

ACDM models predict a faint-end slope for the mass 
function of -1.8 ( [Trentham fc Tully|[2002l ). We measure 
a luminosity function faint-end slope, a ~ —1.28 ± 0.06 
for M81 group galaxies. The steepest faint-end slopes 
consistent with our data are found to be —1.39 ± 0.05 
for our binned differential LF, constructed with all can- 
didates and having presumably overestimated magnitude 
corrections under the assumption of exponential profiles, 
and —1.45 ± 0.06 when fitting a possible faint-end upturn 



observed in the full sample. Even when we apply com- 
pleteness corrections at faint magnitudes, we find faint- 
end slopes flatter than —1.30. When we include only most 
likely member candidates in our LF, we flnd a slightly shal- 
lower slope of —1.27 ±0.06. We therefore expect —1.30 to 
be an upper limit. 

We compare the M81 Group LF with that of other 
nearby groups. Detections of dwarf galaxies in the Local 
Group go much fainter than we recover in M81. However, 
because of the full sky coverage necessary to obtain a com- 
plete census, completeness is still quite low. We therefore 
look only at the LF for Andromeda and its satellites. Large 
areal coverage is still required, including probing into the 
zone of obscuration, but progress is being made. We con- 
struct a cumulative LF for And romeda satellites inclu ding 
all previously known galaxies (van den Bergh"2 006|) and 
recent family additi ons (Ibata et al. 2007: Majewski et al 



20071 



et al. 



Martin et aL]|2006| [Irwin et al . 2008; McGonnachie 



2008) (Figure 35). A fit to the cumulative V band 



counts finds a shallow faint-end slope of — 1.13^Q Qg. This 
is very similar to the be st Schechter function fit for Local 
Group galaxies found by Pritchet & van den Bergh ( 1999 1 
with a — — l.lzbO.l. The recent discoveries of Local Group 
dwarfs have not increased the slope significantly. 

The nearby Cen A group at 3.6 Mpc, in contrast to the 
spiral dominant Local Group and M81 Groups, contains 
a dominant giant elliptical galaxy with radio loud AGN 
and a late type fraction of only ~ 38%. Using B band 
magnitudes for the Cen A group from [Karachentsev et al" 
( j2002b) , we find a faint-end slope a 
a s steep as tha t of M81 (Figure 35). 



1.23r^:^S- nearly 



Boyce et al.[ ( j2001| ) find that steep LFs are generated by 
the combined number counts of dl and dE galaxies, and 
when dis are excluded from the construction of the LF, 
the faint-end slope flattens. This might imply that more 
evolved regions with fewer late types should have shallower 
slopes. Alternatively, from a comparison of poor groups, 
clusters, and field environments, Zabludoff fc Mulchaey[ 



(12000) found that the dwarf-to-giant ratio (and therefore 
the faint end of the LF) increased in regions of higher 
density, and that this increase was specifically due to an 
increase in the fraction of quiescent dwarfs. In our work on 
M81 we find a steeper slope in the dynamically unevolved, 
dl-rich M81 Group than these other two groups and in 
particular the dSph-rich Local Group. However, recent 
work on the luminosity function for two nearby dynami- 
cally evolved groups with large early type galaxy popula- 
tions finds steep faint-end slopes with values of —1.3 ± 0.1 
for the NGC 5846 group (Mahdavi ct al. 2005) and -1.35 
for the NGC 1407 group |Trentham et al.„2006| ), shghtly 
steeper than what we find for M81. These works deter- 
mine group membership on the basis of morphology and 
surface brightness and both reach down to = —12. 

It is possible that the shallower Local Group LF is due 
in part to incompleteness in the surveyed region. In re- 
cent work, Koposov et al. (2007) construct the LF for the 
Local Group including the recent discoveries of ultra-faint 
dwarfs with SDSS. They also include a volume correction 
factor, assuming different number density laws for the dis- 
tribution of halo satellites, and find faint-end slopes rang- 
ing between —1.29 < a < —1.25 for the magnitude range 
— 19 < My < —3, very similar to the slope we find for 
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M81. In this case, with all groups discussed here having 
faint end slopes of a ^ —1.3, we find no clear trend in 
faint end slope as a function of environment within groups 
of galaxies. 

A few recent studies have found steep faint-end slopes 
in group environments. Using a statistical subtraction of 
the background for n earby groups and poor clusters iden- 
tified in SDSS data, Gonzalez et al. ( |2006 ) find a sharp 
upturn in composite LFs at Mr ^ —17 with slopes of 

— 1.9 < a < —1.6. They find this upturn regardless of 
group mass, number of members, or environment. Using 
colors and m orphology to establish membership, |Krusch| 
et al. (20061 also find a steep upturn with a ^ —1.7 start- 
ing around Mb > —15 in Hickson compact groups. While 
we do see a possible increase in slope around Mr' > — 15, 
this appears to be followed by a fiattening at fainter mag- 
nitudes. We do not find any evidence for a slope as steep as 
—1.7. If the slope were this steep in our survey region, we 
would expect to detect the large numbers of bright dwarf 
galaxies implied by this slope. It is hardly conceivable 
that our faint-end slope measurement suffers from selec- 
tion effects compatible with a factor of 10 inconsistency in 
counts at Mr' ^ —10. Down to at least M^/ < —10 we 
do not find such a steep upturn in the poor M81 Group 
when constructing the LF from our nearly complete survey 
of galaxies whose membership is established through more 
direct means. Alternatively, the M81 Group may represent 
a different environment than probed in these two stud- 
ies. Our results are more consistent with the large field 
surveys such as SDSS and COSMOS which find slopes of 

— 1.3 (or —1.5 with corrections for s urface brightness se- 
lection effects) and —1.2 respectively ( Blanton et al.|2005 



Liu et al.|[2008| . With our deeper survey limit we probe a 



fainter dwarf galaxy population yet find a similar slope to 
that found in these redshift studies. 

4.4. Implication for Cosmological Models 

Assuming the LF follows a constant slope of — 1.3, rather 
than turning over at r' > 17.0 where we start to be- 
come incomplete, we would expect to be missing at least 
~ 70 M81 Group members in our survey area in the range 
17<r'<22(-ll< Mr' < -6). If the slope increased 
to —1.8 at faint magnitudes, beginning at A/,' — —12, we 
would expect to be missing over 1700 galaxies in our sur- 
vey area. This is over an order of magnitude greater than 
we would expect from our data. For the intrinsic faint- 
end slope to be this steep, we would have to be missing a 
very large population of very faint and low surface bright- 
ness, low stellar concentration galaxies, and based on our 
simulations, we do not find any evidence that we could 
be missing such large numbers of galaxies brighter than 

Mr' = -7. 

Recently there have been a number of successful at- 
tempts at bringing model predictions closer into line with 
the lower number counts of observations. One idea is 
that the dark matter halo circular velocity is larger than 
extrapolated from measured stellar velocity dispersions. 
This could be caused by dark matter halos having much 
larger extents than the luminous component of the galaxy 
such that velocity dispersions only pr ovide an estimate 
of th e mass in the interior of the halo (Pefiarrubia et al. 
20081. Another explanation posits that dark matter ha- 



los had more mass at higher re dshift but have since lost 
mass through dynamical friction ( Kravtsov et al.|2004 l. In 
both of these scenarios, the observed galaxies today cor- 
respond to larger dark matter halos in the simulations, 
thus alleviating the problem of excess counts predicted at 
the measured luminosities. Another widely used solution 
to the discrepancy incorporates the effects of reionization 
on the smallest mass halos in galaxy formation models. 
Reionization suppresses the accumulation of gas and thus 
star formation, leaving lower mass halos dark. A number 
of recent m odels find good agreement with observations 
dSimon fc Gcha 2007; Kang 2008; Bovill & Ricotti 20q8|. 
jKoposov et aI7[p007[ ) compare the Milky Way LF to 



semi -analytical prediction s and find a reasonable match 



with Benson et al. (2003) who include effects of tidal dis- 
ruption and photoionization with CDM theory predictions. 
Counts to My < —3 agree well within the uncertainties. 
This provides some support for the hypothesis that at least 
part of the cause of the disparity between ACDM predic- 
tions and observations is due to the suppression of gas ac- 
cumulation in small galaxies collapsing post reionization 
although other predictions from their model, such as cen- 
tral surface brightness values, do not agree well with ob- 
servations. In further support of this idea, |Simon fc Geha| 
(20071 compare the observed Local Group circular veloc- 
ity distribution function (a proxy for the mass function) 
to that predicted from the Via Lactea N-body simulation 
(Diemand et al. 20071. They find good agreement when 
assuming that the Local Group halos observed at z=0 were 
the ones which collapsed prior to reionization and which 
correspond to the model objects with the largest values of 
Vcirc at various high redshifts. 

We find a faint-end slope of the M81 Group LF to a 
limiting magnitude of Mr' = — 10 that is much shallower 
than ACDM cosmology would predict but similar to that 
found for Milky Way halo satellites. We do not reach to 
the same ultra faint population as is currently being un- 
covered in the Local Group, and we have not obtained 
kinematical measurements to directly compare the mass 
function with theoretical predictions. The M81 halo has a 
similar mass as the Milky Way, thus solutions to the miss- 
ing halo problem which bring model number counts into 
agreement with observations by associating more massive 
halos with fainter galaxies while lower mass halos remain 
dark could equally well explain the shallow slope we find. 
Therefore, part of the explanation for the 'missing galaxy' 
problem may be in the physics of how mass is converted to 
light. Incorporating the effects of feedback and star forma- 
tion efficiency may help reconcile theory and observations. 
Accounting for the suppression of gas infall into the low 
mass halos of the forming galaxies by reionization in the 
early universe may explain the bulk of the discrepancy and 
naturally explains the noted environmental dependence of 
the LF, with clusters and more dynamically evolved re- 
gions exhibiting slope s with a steeper t han —1.4 (see e.g . 
i Trentham et al.](|2005|) ; [C onsclicc (2002); Chiboucas fc Ma- 
'teo (20061; Ide Propris fc Pritchct (1998); Bernstein et al 



(1995 )) as compared to the shallower slopes typically mea- 
sured in group and field environments, and which we find 
in the M81 Group. 

5. SUMMARY AND CONCLUSIONS 
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We have discovered 22 new candidate dwarf galax- 
ies in the M81 Group from our 65 square degree 
CFHT/MegaCam survey, designed to extend beyond the 
second turnaround radius for the group. Of these candi- 
dates, we beHeve 12 are Hkely members. The remainder 
may prove to be background galaxies, foreground associa- 
tions of stars, or artifacts in our data. The 12 likely mem- 
bers consist of 3 candidate BCDs lying near the periphery 
of the group, 1 tidal dwarf candidate lying within the HI 
tidal bridge between M81 and NGC 3077, and a number 
of dl and dSph candidates. No new large diffuse or any 
ultra compact dwarfs are detected. The large fraction of 
late to early types previously noted in the M81 Group is 
reconfirmed with our new objects. Assuming all objects lie 
at the distance of M81, we find r' absolute magnitudes for 
these new objects range from -6.7 to -13.3. The faintest 
object, if real, has a measured size of 90 pc that would 
cause it to encroach on the 40 - 100 pc size gap region 
between globular clusters and dwarf galaxies. From false 
star /galaxy tests we expect to detect nearly 100% of the 
group member dwarfs down to Mr' = — 10 and over 50% 
to Mr' = —9, not quite into the regime of the recently 
discovered ultra faint Local Group galaxies. Including all 
22 previously known M81 Group members in our survey 
region, we construct the group differential and cumula- 
tive luminosity functions and find modestly steep faint-end 
slopes of ~ —1.30 ± 0.06. Including only the most likely 
members, we find a similar slope of —1.27 ± 0.06. This 



slope is steeper than what has been found in the case of 
the Andromeda satellites, but may be consistent with the 
Local Group after taking into account all new discover- 
ies and including corrections for completeness. Even with 
the addition of the 22 new candidate dwarf galaxies in the 
M81 Group, number counts remain an. order of magnitude 
below cosmological predictions of halo counts. 
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^Bt are estimated by eye on POSS-II 

^angular dimensions correspond to an isophotc 28m arcsec 



Table 2 

Spectroscopy of 5 dwarfs in M81 
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Table 4 



Photometry of previously known M81 group members 



Name 


Q 


5 




typo 


e 


Mo 
T arcscc 


r ai cscc 


Re 
arcsec 




K' 
mag 


r' ^ 


'■cor 

mag 


M,,, 

cors 


Hoi 


9 40 28.2 


+71 11 


11 


dl 




25.0 


25.2 


66.7 


13.9 


0.1 


13.4 


13.4 


-14.4 


F8D1 


9 44 50.0 


+67 28 


32 


dSph 




26.5 


26.6 


67.3 


15.3 


0.2 


14.0 


14.7 


-13.0 


FMl 


9 45 10.0 


+68 45 


54 


dSph 


0.14 


24.8 


25.2 


21.1 


16.5 


0.2 


15.9 


16.2 


-11.5 


N2976 


9 47 15.6 


+67 54 


49 


Sc pec 


0.60 


20.4 


20.7 


43.6 


10.5 


0.2 


9.7 


9.8 


-18.0 


KK77 


9 50 10.0 


+67 30 


24 


dSph 


0.29 


24.5 


25.2 


36.7 


15.4 


0.4 


14.7 


14.9 


-12.8 


BK3N 


9 53 48.5 


+68 58 


09 


dl 


0.01 


25.9 


24.8 


6.7 


18.6 


0.2 


17.8 


17.8 


-10.2 


M81 


9 55 33.5 


+69 04 


00 


Sb 




17.1 


19.3 


81.1 


7.8 


0.2 


7.2 


7.2 


-20.6 


M82 


9 55 53.9 


+69 40 


57 


In- 




17.4 


19.1 


60.0 


8.4 


0.4 


7.9 


7.9 


-19.9 


K61 


9 57 02.7 


+68 35 


30 


dSph 


0.25 


23.9 


24.7 


38.4 


14.8 


0.2 


14.5 


14.5 


-13.3 


A0952+69 


9 57 29.0 


+69 16 


20 


pec 










16.3 


0.2 








HoIX 


9 57 32.4 


+69 02 


35 


dl 


0.24 


24.1 


24.7 


46.8 


14.3 


0.2 


13.9 


14.0 


-13.8 


N3077 


10 03 21.0 


+68 44 


02 


Irr 


0.24 


19.6 


20.2 


37.5 


10.3 


0.2 


10.0 


10.1 


-17.8 


Garland 


10 03 42.0 


+68 41 


36 


pec 










16.5 


0.2 








BK5N 


10 04 40.3 


+68 15 


20 


dSph 


0.20 


24.5 


25.1 


22.1 


16.4 


0.1 


16.1 


16.1 


-11.8 


K63 


10 05 07.3 


+66 33 


18 


dSph 


0.18 


23.9 


24.1 


24.9 


15.1 


0.2 


14.4 


14.8 


-13.0 


K64 


10 07 01.9 


+67 49 


39 


dSph 


0.33 


23.2 


23.8 


22.1 


15.1 


0.2 


14.6 


14.7 


-13.1 


IKN 


10 08 05.9 


+68 23 


57 


dSph 










(15.7) 


0.2 


(15.5) 


(15.5) 


(-12.3) 


HS117 


10 21 25.2 


+71 06 


58 


dl 


0.39 


24.3 


24.8 


19.9 


16.3 


0.3 


15.5 


15.9 


-11.9 


DD078 


10 26 27.9 


+67 39 


24 


dSph 


(0.10) 


(24.8) 


(25.0) 


(25.0) 


(15.9) 


0.1 


(15.0) 


(15.0) 


(-12.8) 


IC2574 


10 28 22.4 


+68 24 


58 


SABm 




22.6 


24.0 


142.0 


11.5 


0.3 


10.2 


10.3 


-17.7 


DD082 


10 30 35.0 


+70 37 


10 


Im 


0.40 


21.8 


22.7 


33.9 


13.1 


0.1 


12.9 


12.9 


-15.1 


BK6N 


10 34 31.9 


+66 00 


42 


dSph 


0.32 


24.6 


24.6 


16.6 


16.2 


0.0 


15.9 


15.9 


-12.0 



Note. — Values in parentheses highly uneertain due to presence of bright foreground stars. Values for brightest galaxies may be underestimates due 
to sky subtraction errors in the TERAPIX reduction pipeline. 

^Corrected for extinction and flux lost in sky assuming an exponential profile. 

"^Corrected for extinction and flux lost in sky assuming a Sersic profile. 

'^Absolute Sersic profile corrected magnitude, using distance moduli from Table 2 of |Karachentsev et a"L] l |2002a[ l or assuming 27.8 where not provided. 



Table 5 



Faint-end slope from differential luminosity functions 



Sample 




Q (^^TTieas) 


" «or(.xp 


)) 


" ('■coT-(ser=ic 


)) 


All (completeness 


corr.) 


-1 29+°"* 


-1.28t°- 


04 
.05 


-1.32+°- 


04 
06 


All 




-1 26+°"" 


-1.39+0 


05 
.05 


-1.36+^- 


09 
OS 


Best 




-1 20+°-^" 


-1.34+«- 


07 
.07 


-1.30+;;- 


10 
10 


Original 




-1 19+0"^ 


-1.24i«„ 


09 

,08 


-1.26+°- 


08 
08 


An (-15 < r' < - 


-10) 


-1 43+0"3 


-1.44t«„ 


07 
,07 


-1.46i°: 


10 
10 



Table 6 



Faint-end slope from integrated luminosity functions 



Sample 


OL (r^ 




" «or(.xp)) 


" «or(scrsic 


)) 


All (completeness corr.) 


-1.29 


+ 0,02 
-0,05 


1 00 + 0,03 
-l-^8„o 05 


-1.28+° 


02 
,05 


All 


-1.29 


+ 0,04 
-0,06 


, r,Q + 0,05 

-^■■'''-0.08 


-1.28+0' 


04 
.04 


Best 


-1.27 


+ 0,05 
-0,06 


, ,7+0,04 
-^■^'-0-07 


-1.26+;5' 


05 
,04 


Original 


-1.26 


+ 0,07 
-0,10 


-1 24+"'"' 
1.^^ 0.07 


-1.25i° 


06 
,06 
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40 20 Q 

Supergaloctic Longitude 

Fig. 1. — The extended region around M81. All galaxies with velocities in the Galactic standard of rest less than 400 km s^^ in an area 
of 60 X 40 degrees in supergalactic coordinates. Circles: spirals and irregulars. Squares: early types. Large, mid, and small symbols identify 
galaxies with Mg < —20, —20 < Mg < —16, and A{g > —16, respectively. The box outlines the area of the CFHT M egaCam survey and 
the small open triangles identify the newly discovered dwarfs. Contours illustrate the dust maps of jSchlegel et alT| | |1998| at intervals of 0.2 r' 
magnitudes of extinction. 
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Fig. 2. — Histogram displaying the seeing range for all combined images in our survey. 
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Fig. 3. — Tests of photometric accuracy. We compare our r' total magnitudes to transformed R magnitudes of [Perelmuter fc Racine] l |1995[ |. 
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Fig. 4. — Thumbnails of the 3 new M81 candidate BCDs. Images are 1.5 arcmin on a side. 
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Fig. 5. — We display a mosaic of images of 9 new candidates that were detected by eye and/or 2-point correlation of resolved points. Images 
are 1.5 arcmin on a side. 
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Fig. 6. — While these six candidates have a slight degree of resolution, many of these may prove to lie in the background of the M81 Group. 
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Fig. 7. — Four candidates which may turn out to be artifacts (excess noise or foreground concentrations of stars) or distant galaxy clusters. 
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Fig. 8. — Map of our M81 survey region (large box almost filling the figure). Circles indicate previously known group members while stars 
denote the location of the new group candidates. The large dashed circle is the projection of the putative surface of second turnaround at 
230 kpc from M81. 
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Fk;. 9. Four examples of object concentrations detected with our 2-point correlation routine that we believe to be distant galaxy clusters. 
Image sections are 2.2 arcmin across. 
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r 



Fig. 10. — Wc display histograms of the resolved objects in the new M81 candidates. Number counts of objects recovered by SExtractor 
within a 30-45 arcscc radius (sec text) centered on the new candidates are binned as a function of magnitude. The top row includes 4 objects, 
including 3 BCD candidates, with the largest resolved stellar populations. The next 2 rows include a further 8 good candidates. The 4th row 
contains 4 objects that may prove to be artifacts, while the bottom 2 rows include possible more distant galaxies. The shaded histograms 
represent number counts in an identical sized region 2.5 arcmin west of each galaxy. 
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Fig. 11. — Examples of artificial galaxies that were added to real images with which we test the recovery efficiency of our detection 
techniques. Images are 1.5 arcmin on a side. 



Discovery of New Dwarf Galaxies in M81 



29 



0.5 (kpc) 1 1.5 




Fig. 12. — Results of 1200 artificial galaxy recovery tests. Small pentagons represent all artificial galaxies. Open pentagons are those 
galaxies that landed at least partially in chip gap regions of the MegaCam field or are largely obscured by bright saturated stars or other 
objects and are less likely to be recovered. Stars around these points denote that the simulated galaxy was recovered. Overlaid as black 
diamonds are the candidate real galaxies in this magnitude - effective radius plane. The solid line denotes the magnitude at which the recovery 
drops to 50% (in 2-magnitude bins), while the dashed line marks the boundary of 90% completeness. 
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Fig. 13. — We plot the recovery of 1200 artificial galaxies in the magnitude-surface brightness plane. As before, pentagons represent all 
artificial galaxies. Stars around these points indicate that the simulated galaxy was recovered. Overlaid as black (gray) diamonds are the 
candidate real galaxies (BCDs) in this magnitude - effective surface brightness plane. Circles represent previously known M81 Group galaxies. 
DDO 82, shaded in gray, is a potential BCD. Lines of constant effective radius for exponential profiles arc labeled. We draw a solid curved 
line which bounds the region of recovered artificial galaxies. 
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Fig. 14. — Recovery results for objects simulated with 0.5 magnitudes of extinction. Symbols as in previous figures. The solid line represents 
the magnitude at which the recovery drops to 50% for these galaxies, while the dashed line denotes the magnitude at which the recovery for 
artificial galaxies with 0.-0.2 magnitudes of extinction drops to 50%. 
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Fig. 15. — Recovery of 400 galaxies added with 0.9 arcsec seeing. Symbols as in previous figures. We compare results to that of the original 
1200 simulations. The dashed line represents the magnitude at which the completeness drops to 50% for artificial galaxies added with seeing 
less that 0.75 arcsec. The solid line displays the 50% completeness for artificial galaxies added with 0.9 arcsec seeing. 




Fig. 16. — Spectra of 4 M81 group members obtained with Subciru/FOCAS. All spectra are shifted to the rest frame given the cross- 
correlation measured velocities as labeled. 
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d0958 + 66 6-m telescope 2006/20/11 
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Fig. 18. — Spectrum of d0958+66 from the 6m BTA telescope. Strong Ha emission is detected and used to measure a heliocentric radial 
velocity. 
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Fig. 19.— 
S-knots. 



Spectrum of dl028+70 from the 6m BTA telescope. Ha emission is detected in two separate knots which we label as N and 
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Fig. 20.— 



Color images based on 2 color photometry of 9 dwarf candidates. 
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Fig. 21. — Curve of growth (top) and surface brightness profiles (bottom) for our 4 brightest candidates with best cumulative (top) and 
surface brightness profile (bottom) Sersic function fits overlaid. In the bottom plots, we also overlay the best cumulative Sersic function fit 
(green long dash) and best exponential surface brightness profile fit to outer radii only (cyan short dash). 
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Fig. 22. — Fits for 4 of our candidate galaxies, as in Figure [21] Curve of growth profile data points which exhibit a sharp increase in slope 
beyond 10 arcsec were truncated during profile fitting. For d0926+70 and d0934+70 this rise was determined to be due to nearby bright 
cirrus. For d0959+68 we were unable to obtain a surface brightness profile fit due to the fact that this object consists almost entirely of bright 
stars without a significant low surface brightness component. 
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Fig. 23. — Fits for 4 of our candidate galaxies, as in Figure [211 
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Fig. 24. — Fits for 3 potential background galaxies, as in Figure [21] 
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Fig. 25. — Fits for a further 3 likely background objects, as in Figure [21] 
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Fig. 26. — Fits for 4 of our poorer candidates, as in Figure [2T] 
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Fig. 27. — Differential luminosity function for the M81 Group, using measured magnitudes (top), total magnitudes corrected assuming 
exponential profiles (middle), and total magnitudes assuming Sersic profiles (bottom). Counts are normalized to the area of the survey 
coverage. Solid red circles include previously known group members and all candidates from this survey. Open circles are corrected for 
incompleteness. Green triangles include only the best candidates from this work. Small blue circles represent only previously known group 
members. The hatched region denotes where this survey is incomplete; brightward of this we expect nearly 100% completeness in the detection 
of member dwarfs. The best power law fits to the data brightward of our completeness limit are shown by the dashed lines, while the solid 
line is the best fit for corrected counts. Dotted lines are faint-end only fits to all counts within the magnitude range — 15 < r' < —10. The 
arrow in the final data point signifies that a minimum completeness correction based on good seeing/low extinction simulations was applied. 
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Fig. 28. — We plot detection efficiency as a function of magnitude. The solid line is the best fit to the good seeing, low extinction case, 
while the dashed line incorporates the survey average effects of seeing, extinction, single exposure regions, and loss of survey area due to 
coverage by brighter objects. 
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Fig. 29. — Cumulative luminosity function for the M81 Group, using measured magnitudes (top), magnitudes corrected assuming exponential 
profiles (middle), and corrected magnitudes assuming Sersic profiles (bottom). Counts are normalized to the area of the survey coverage. Large 
solid circles include previously known group members and all candidates from this survey. Open circles are corrected for incompleteness. The 
hatched region denotes where this survey is incomplete; brightward of this we expect nearly 100% completeness in the detection of member 
dwarfs. The best cumulative Schechter function fits to the data brightward of our completeness limit are shown by the dashed lines. The 
solid red line is the best fit to completeness corrected counts. Dotted lines are zblcr fits for the sample of 'best' candidates. 
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Fig. 30. — We display 1 u error ellipses for the luminosity function parameters a and M* for our best (solid) and all (dashed) candidate 
samples using Sersic corrected magnitudes. The values for M* are unconstrained at the bright end. 
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Fig. 31. — We display the M81 HI map of |Yun et al.| l [l994[ |. Overlaid (white dots) are the positions of the five suspected tidal dwarfs. 
d0959+68 is found to lie in the tidal bridge between M81 and NGC3077. Adapted by permission of Macmillan Publishers Ltd: Nature, 
copyright 1994. 
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Fig. 32. — Locations of our M81 candidates on a total magnitude vs half-light r adius plane (filled circles), along with previously known M81 
alaxies (filled squares), Local Group previously known galaxies (open squares) Jlrwin fc Hatzidirnitriou|1995l rMateo | 1998| | McConn^ 



Irw in 2006 ) and recent discove ries (open circ les 
leTal. 2006 Martin ct al. 2006| |Ibata et al!|20 07 



Simon & Geha 2007 Belokurov ct al. 2007a ;'WiTIman ct al 



. 20^ 



[Walsh et al.|2DD7 



IZucker 



Irwin ct al. 2008). 'i'o transform V magnitudes to r', we assume'F' = V — U.84(V — H 
'iikugita ct al. f996b). i'he rough locations of globular clusters and Fornax UCDs on this plane are shown. Globular Clusters have sizes 
at range up to ~ 40 pc, while Local Group galaxies have sizes generally larger than ^ 100 pc. The size gap region falls between the two 
dashed lines. Lines of constant effective surface brightness are also shown. 
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Fig. 33. — We display the M81 galaxy detections superimposed on the |SchlegeI et al.| | |1998| dust map. Red diamonds represent newly 
discovered candidate group members while yellow crosses correspond to previously Jtnown group members. 
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Fig. 34. — Projected distribution of early type (red hexagons) and late type (blue circles) galaxies. Green squares indicate candidate tidal 
dwarfs. The location of BCD candidates, including DDO 82, are marked by blue rings. We include in this plot only previously known and 
best candidates. Size indicates surface brightness with larger points denoting brighter effective surface brightness. 
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Fig. 35. — Cumulative LFs for Cen A {B band, squares), Andromeda {V band, stars), and M81 (r' band, open circles) satellites. For clarity, 
we do not include error bars. Applying a cumulative Schechter function to these data, we find faint-end slopes of — l.lS^Q Qg for Andromeda 
systems, -1.23±tl'°t ^en A, and -1.28lo;og 



